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We discuss latest results of lattice QCD simulations with dynamical fermions. Special emphasis is paid to the 
subjects of the static quark potential, the light hadron spectrum, T spectrum, and the pion-nucleon-sigma term. 
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1. INTRODUCTION 

The investigation of the influence of sea quarks 
on the physics of strong interaction has been a 
field of vital interest ever since the formulation of 
lattice quantum chromo dynamics. 

However, the computer simulation of full QCD 
processes is a very challenging task, and it is 
only in these days, that computers and algorithms 
might become powerful enough to produce results 
which can be connected reliably to continuum 
physics. 

Several large scale lattice simulations of full 
QCD have been started recently, unfortunately 
most of them being still in the phase of lattice 
calibration or program development. The CP- 
PACS, UKQCD and QCDSP projects belong to 
this category. Status reports and first (prelimi- 
nary) results can be found in |],|, §,[§, §]. The 
bulk of - at least semi-final - results achieved 
this year comes from the SESAM and the T%L 
|6| collaborations. Therefore, these results will 
serve as a guideline for this talk. In the first sec- 
tion I will discuss the status of static potential 
measurements. The second section deals with 
slowly moving quarks, i.e. the NRQCD results 
come into focus. The third section summarizes 
the SESAM/T%L investigations on the light spec- 
trum, and finally the lattice measurement of the 
pion-nucleon-sigma term ajv is reviewed in sec- 
tion four. R 



quenched QCD, is that the energy string between 
a quark and it's anti-quark should break once the 
distance between these particles becomes large 
enough to produce a quark anti-quark pair out 
of the vacuum. Clearly, one would like to verify 
this important quality in full QCD lattice simu- 
lations. One prominent quantity, which should 
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Figure 1. SESAM/T X L collaboration. The 
rescaled static potential in full and in quenched 
QCD. The data at k = 0.1575 comes from ana- 
lysis of the T\L lattice. The methods of analysis 
have been explained in detail in |^/ . 



2. STATIC POTENTIAL 

One of the most exciting features that we ex- 
pect to be a property of full QCD, but not of 

1 For flavour dependence of chiral quantities sec|10|]l l|. 



exhibit a clear sign of string breaking, is the po- 
tential between static quarks. The signature that 
we expect is a flattening out of the potential at 
a given distance rv ea fc. Unfortunately, this has 
not been observed yet in lattice QCD simulations 
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up to distances of r ~ l.bfm and quark masses 
down to approximately the strange quark mass. 
The SESAM/T^L collaboration has performed a 
high precision measurement of the static potential 
fl2|| with n/ = 2 dynamical Wilson fermions at 
(3 = 5.6. SESAM has used a lattice size 16 3 x 32 at 
four different values of the sea quark mass, n sea = 
0.1560,0.1565,0.1570,0.1575, corresponding to 
m^/mp = 0.839(4), 0.807|,0.755(7),0.69(1), with 
200 statistically independent gauge configura- 
tions. The T^L collaboration has used lattice 
size 24 3 x 40, two sea quark masses, K sea = 
0.1575,0.1580, corresponding to ra^/rap ~ 
0.69,0.55, and a statistics of presently ~ 80 con- 
figurations. The result - at two different values 
of the quark mass - is shown in fig.[l]. Obviously, 
there is still no clear sign of string breaking in 
this range of quark masses and up to a physical 
distance of r ~ 2.0/to. 
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Figure 2. UKQCD collaboration (prelimi- 

nary). The rescaled static potential in full and in 
quenched (/3 = 5.7) QCD. r/rn = 3 corresponds 
roughly to r ~ 1.5 fm. The naive expectation of 
string breaking is displayed at the right edge. For 
details see J^/ . 



The result of SESAM/T^L is supported by 
the (preliminary) findings of the UKQCD collab- 

2 Configurations with this sea quark mass are currently 
analyzed. 



oration j3J]. Their result, obtained with an im- 
proved clover action [|l3| at rif = 2, (3 = 5.2, 
and a lattice volume of 12 3 x 24 is displayed 
in fig.^. In order to guess, at which energy 
V(rbreak) the string should break, UKQCD has 
also calculated the mass of the static light meson 
(Msl) at Kught = 0.1390. Naively, one expects 
V(rb rea k) — 2 x Msl- The data in fig.[| however, 
doesn't seem to obey this naive expectation, al- 
though there is still a small chance that they may 
do so (at K sea = 0.1390). This uncertainty can be 
removed soon, as UKQCD's simulation proceeds. 

So, finally one has to tackle the question why 
string breaking could not be observed yet in full 
QCD determinations of the static potential. One 
explanation could be that currently available sea 
quarks are still too heavy. The energy necessary 
to produce a quark anti-quark pair out of the 
vacuum would then be available only at larger 
distances. This is however in sharp contrast to 
the naive expectation and would therefore require 
further investigation. 

Another explanation, which appears much 
more likely to the author, is that the Wilson loop 
operator, which is used to calculate the potential, 
does not have sufficient overlap with the state in 
which the string is broken. In order to corrob- 
orate this scenario one would have to compare 
the potential energy with the ground state energy 
of an operator, which is designed to describe the 
physical situation of a broken string. Such an op- 
erator could be that of a system of two static- light 
mesons. Its ground state would have to be eval- 
uated with respect to the spatial distance of the 
mesons. The crossover point between V(r) and 
the ground state energy of this operator would 
then determine r^reak- 

Before we turn to the discussion of NRQCD re- 
sults, we would like to emphasize that the influ- 
ence of dynamical fermions can be seen clearly in 
the short distance part of the static potential, i.e. 
in its coulombic contribution. Fig. |^ compares 
the values of the coulomb coefficient e at several 
sea quark masses and in the quenched approx- 
imation. The difference between quenched and 
full QCD results is statistically significant and re- 
veals that the coulombic part becomes stronger in 
full QCD. 
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Figure 3. SESAM/TxL collaboration. The 
coulomb coefficient e in full and in quenched ((3 = 
6.0J QCD. The analysis has been done analogous 
to§. 



3. T Spectrum 

The determination of heavy particle proper- 
ties with lattice methods is a notoriously diffi- 
cult task, as large cutoff effects may contaminate 
the results. A very promising attempt to get rid 
of these contaminations is the lattice formulation 
of NRQCD, which, in principle, is able to con- 
trol its systematic errors. As the computational 
effort needed to exploit NRQCD is small com- 
pared to standard lattice simulations, it appears 
to be an ideal laboratory to study the influence 
of sea quarks on the spectrum of heavy mesons 
with high statistical precision. 

It has been reported last year that the the 
T spectrum energy levels are shifted towards their 
experimental values once sea quarks are switched 
on. This behaviour had been seen in the spin 
independent as well as in the spin dependent parts 
of the spectrum. 

Encouraged by this result the SESAM collabo- 
ration has now performed an NRQCD analysis of 
the T spectrum with n/ = 2 Wilson fermions|Q. 
With this data in hand, which we display in fig. 
[|, the situation now looks much less definite. Al- 
though SESAM consolidates the last years pic- 
ture with respect to the spin independent part of 
the spectrum, there is no significant effect of sea 
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Figure 4. Spin independent (upper graph) and 
spin dependent parts of the T spectrum. Open 
circles: NRQCD collab. quenched^; filled cir- 
cles: NRQCD collab, nj = 2, KS fermions; stars: 
SESAM, quenched^, crosses: SESAM, n f = 2. 
Wilson fermions, M: Manke, nf — 0/uT 



quarks in the spin dependent part. 

In order to explain the difference between the 
results of the NRQCD collaboration |J] and the 
SESAM results, one has to remember that the 
definition of the NRQCD Hamiltonian is by far 
not unique. The NRQCD collaboration has in- 
cluded only terms up to 0(Mv 4 ) and used the 
plaquette definition of the average link for tad- 
pole improvement. SESAM went up to 0(Mv 6 ) 
and used the Landau gauge definition of the av- 
erage link for tadpole improvement. The effects 
of such changes have been studied in detail by 
]f^ . p^ | within the quenched approximation. For 
the spin dependent splittings of the T system they 
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find them to be of the order 10% - 20%. We con- 
clude that, using state of the art NRQCD, un- 
quenching effects in the spin dependent part may 
be hidden by systematic uncertainties. 

4. LIGHT SPECTRUM 

The accurate calculation of the light hadron 
spectrum, including the strange quark sector, is 
still an open problem for full lattice QCD. Apart 
from the lack of statistics, finite size and finite 
cutoff effects, one reason is that one cannot in- 
clude light and strange sea quarks at the same 
time yet. 

Thus, the usual way of analyzing spectrum 
data has been to look at light hadrons in a sea 
of light quarks only, whereas the properties of 
hadrons containing strange quarks have been ap- 
proximated by use of a sea of strange quarks. 

The SESAM collaboration has now developed 
a more general type of analysis, which goes be- 
yond the restriction of equal values of sea and 
valence quark mass. With that method one ex- 
ploits a given full QCD gauge configuration - 
with a definite value of the sea quark mass m sea 
- with respect to several valence quark masses 
m va i. Guided by chiral perturbation theory one 
can then parameterize hadron mass results with 
respect to both m sea and m va i. A simultaneous 
fit to all hadron data finally allows to access arbi- 
trary points in the (m sea ,m va i) plane. Of course, 
the validity of this ansatz has to be checked care- 
fully with the data. 

SESAM has applied this method to the set of 3 
(sea quarks) x 5 (valence quarks) combinations, 
with a high statistics of 200 independent gauge 
configurations per sea quark, using a linear ansatz 
both in m va i and m sea . The valence quark masses 
are located in the same range as the sea quark 
masses (see above) . Fig. || illustrates the quality 
of the approach for the case of pseudo-scalar and 
vector mesons. Although the data is reasonably 
well described with a linear ansatz, one cannot 
rule out quadratic terms in m q at this stage of 
the analysis. In order to check on this SESAM is 
currently collecting data at one additional value 

3 We thank M. Liischer for pushing our thoughts into this 
direction. 



of the sea quark mass. 

With the current setup - 3 values of dynami- 
cal quarks and linear parameterization - SESAM 
finds no significant unquenching effect in the light 
hadron spectrum. This is similar to the findings 
of the MILC collaboration jl(| , which uses stag- 
gered sea quarks. 

Finite size effects are a potentially dangerous 
source of contamination, especially if one tries 
to push lattice masses towards the chiral point. 
As almost nothing is known about the size of 
these effects in the case of dynamical Wilson 
fermions|l8|], the T%L collaboration has extended 
the SESAM lattice to size 24 3 x 40. The sim- 
ulation is now being performed at the smallest 
SESAM sea quark mass , k = 0.1575, and at an 
even smaller mass, k = 0.1580. With the current 
statistics - 80(200) TxL(SESAM) configurations 
- a direct comparison reveals finite size effects of 
0(5%), but at this stage of the simulation it is 
clearly too early to draw definite conclusions. 

5. PION-NUCLEON-SIGMA TERM 

The determination of the pion-nucleon-sigma 
term 

a N = mt{P\uu + dd\P) , m/ = i (m u + m d ) (1) 

has a long and eventful history, reaching back 
into the late seventies Although experimentally 
as well as theoretically difficult to determine, it 
is an important quantity, measuring chiral sym- 
metry breaking in QCD. It is also a fascinating 
quantity, as it is expected that pure vacuum fluc- 
tuations make up at least half of its size. Since 
we are looking for sea quark effects, it appears 
mandatory to study in full lattice QCD. 

The cheap way to calculate erjv in full lattice 
QCD is simply to determine the proton mass at 
several quark masses and to apply the Feynman- 
Hellman theorem 

mi(P\uu + dd\P) = mj^-^ . (2) 
ami 

However, in a statistical simulation, the outcome 
of this method may be rather unprecise: suppose 
one would be able to determine mp(m q ) to a pre- 
cision which allows to resolve the quadratic term 
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Figure 5. SESAM/T\L collaboration. Meson masses as a function of quark masses. m & q >' — 
1/2(1/ K e ff ~~ 1/ K c); K c — 0.15846(5). Trig'* is the average valence quark mass of the meson. Sym- 
bols: (ss) data: * and solid lines; (sv) data: O and dashed-dottet lines ; (vv) data: □ and dashed lines 



of this function. Then, because of eq.g, the quark 
mass dependence of a n would be known only to 
the order linear in m q . Thus, in order to extract 
maximal information from a given set of configu- 
rations, it seems advantageous to choose the hard 
way, i.e. to calculate the connected and the dis- 
connected parts of eqjlj 

The SESAM collaboration has followed this 
way of analysis. The connected part has been 
determined with the standard procedures. The 
disconnected contribution to o~n is given by the 
correlation between the proton propagator and 
the quark loops. The latter were calculated by ap- 
plication of the stochastic estimator method with 
complex Zi noise. The statistical quality of the 
correlation signals could be increased substan- 
tially by using a plateau-like sampling of quark 
loops |H| instead of the usual summation over all 
quark loops in space and time. 

Fig.|| displays the SESAM result for the (un- 
renormalized) amplitude (P\uu + dd\P). Both 
contributions, connected and disconnected, in- 
crease as the quark mass decreases. A linear ex- 



trapolation to the light quark mass yields 

o N = 19.65(4.39)MeV . (3) 

Here we have multiplied with the light quark 
mass, mi = 0.5(1/ ki — 1/k c ), ki = 0.158415, k c = 
0.158458, and with the inverse lattice spacing, 
a -1 = 2.33GeV. This value (eq.||) is almost a fac- 
tor two smaller than the results obtained from 
phenomenology and from quenched QCDpH 
p2[ . We emphasize that no renormalization factor 
is involved here, as cat is a renormalization group 
invariant quantity. 

In order to find out where this large deviation 
comes from, we compare in fig. @ the SESAM re- 
sults with previous quenched results. Obviously 
all amplitudes are in reasonable agreement. Thus 
we conclude that the large deviation in o~n is 
caused by the large difference between quenched 
and unquenched quark masses^], found in |^4|,|l^] . 
It is an open issue pq|, whether present lattice 
techniques are sufficient to determine mi reliably 
in full QCD. But lattice calculations of cat appear 



4 This has been anticipated by M. Okawa| 
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Figure 6. SESAM collaboration. Connected and 
disconnected amplitude (P\uu + dd\P) as a func- 
tion of the quark mass. The extrapolated values 
are depicted at m q = . 
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Figure 7. Connected (left) and disconnected con- 
tributions to (P\uu + dd\P) . The data is taken 
from fli|jgi[ / and has been renormalized with the 
standard tadpole improved renormalization factor 
Z s ■ 



to exclude already that uncertainties due to mass 
renormalization factors could be responsible. 

6. CONCLUSIONS 

Sea quarks are shy objects. The only places 
where we could at least see the footprints of Wil- 
son fermions are the short distance part of the 
static potential and the spin independent part 
of the T spectrum. Obviously very accurate 
simulations, very refined analysis methods and 
a clear understanding of the underlying physics 
are needed to observe sea quarks at work. The 
first promising steps in this direction have been 
done recently, and we are looking forward to the 
next generation of full QCD simulations, namely 
UKQCD, CP-PACS and QCDSP. 
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